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(p ■ Abstract 



C^ . We have developed a class of idealized models of chromonic molecules which 

are miscible in water, but which can form aggregates which in turn organize 
C^ ■ 

Vh ■ into lyotropic liquid crystal (LC) phases. By carrying out Monte Carlo simu- 

P^ . lation in a binary mixture of model chromonic and water molecules, we have 

o, 

, ^,. studied the effect of concentration and molecular shape on the nature of re- 
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suiting mesophases. We have also computed the free energy associated with 
the formation of chromonic columnar aggregates by umbrella sampling. Our 



> 
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f~^ . free energy computation helps us to verify the isodesmic behavior which is 

characteristic of chromonic systems. For isodesmic behavior the addition of 

j^ . each chromonic molecule to the columnar aggregates is accompanied by a con- 

stant free energy increment and the net free energy decreases with increasing 

G , column length. 
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I. INTRODUCTION 

Recently there is a growing interest in a new class of lyotropic liquid crystals (LCs), so 
called lyotropic chromonic liquid crystals (LCLCs) due to their diverse potential application 
in many areas such as dyes, drugs, nucleic acids, antibiotics and anti-cancer agents [0],. 
Dichroic thin films formed from LCLCs can be used both as alignment layers and polarizers 



for LC cells 0. These materials are fundamentally different from conventional amphiphilic 
systems: their molecules are disc-like or plank-like rather than rod-like, rigid rather than 
flexible and aromatic rather than aliphatic. Due to their disc-like molecular shape chromonic 
molecules stack face to face forming columnar aggregates. Though there is no clear under- 
standing of the nature of forces responsible for this face-to-face stacking, it is generally 
believed that the vr-vr interaction of the aromatic cores is the main mechanism of molecular 
face-to-face stacking [^ . Hydrophilic ionic groups at the periphery of the molecules make 
them water soluble. The two principal chromonic mesophases are chromonic A^ and M 
phase. The basic structural unit of both phases is the untilted stack of the molecules. The 
A^ phase is the nematic phase in which the molecular stacks has orientational order but no 
positional ordering and is formed when the concentration is low. At higher concentration 
these molecular stacks arrange themselves in a hexagonal pattern forming the M phase. 
In recent years, many new chromonic phases has been reported 0,^ which are formed by 
aqueous cyanine dyes. Some of the liquid-crystalline aggregates formed by the cyanine dyes 
are shown in Figure [I|. 

Unlike the conventional amphiphilic system, the chromonic system does not have any 
critical micellar concentration (CMC). To avoid the unfavorable contact of the hydropho- 
bic part of the molecules with the solvent, the molecules start to stack together forming 
columns even in dilute solution. As the column length increases, the fraction of the total 
molecular surface exposed to the solvent decreases. There is widespread experimental evi- 
dence that there is no optimum column length or aggregation number like the conventional 
amphiphilic system forming micellar aggregates. The micelle has a optimum size depending 
on the packing parameters which in turn depend on the ratio of the size of head and tail 
group 0. The chromonic aggregation behavior is termed as 'isodesmic' in contrast to the 
'non-isodesmic' behavior observed for simple amphiphilic system. For isodesmic behavior, 
addition of each chromonic molecule to the columnar aggregate is accompanied by a constant 
free energy increment. Therefore, with the increase in column length the free energy grad- 
ually decreases, in contrast with the non-isodesmic behavior of simple surfactant system. 



where micelle formation is accompanied by sharp free energy minima. The self-assembly of 
chromonic molecules is governed by many factors such as molecular structure, concentration, 
temperature, solvent polarity and ionic strength. So far there is no clear understanding of 
the relationship between molecular structure and the structure of the supramolecular ag- 
gregates formed by the LCLC materials. Consequently many different aggregate structures 
have been proposed based on different experimental results. 

In this paper, we have developed idealized molecular models of LCLCs which capture 
the overall shape of the LCLC molecules and the nature of their interactions in an approx- 
imate way. To understand the self-assembly of chromonic molecules, we have carried out 
simulations of binary mixture of LCLC molecules and water, and have studied the effect of 
concentration, temperature and molecular shape on the nature of the resulting mesophases. 
We have also computed the free energy associated with the addition of each chromonic 
molecule to a columnar aggregate. The method gives us the possibility to test the widely 
accepted isodesmic behavior of the chromonic aggregation. 

The paper is organized as follows: in section 2 we briefly describe our model and the 
simulation technique. The results on the columnar aggregates, in bulk water, formed by the 
chromonic molecules are reported in section 3. We describe the method of computation of 
free-energy by umbrella sampling in section 4. Finally, a summary of the main results and 
the conclusions drawn from these are given in section 5. 

II. THE MODEL 

We have modeled chromonic molecules as diamond shaped constituted of 9 (model 1) 
tangent spheres of diameter a bonded together (figure 2). We also considered another 
variant of the model where the chromonic molecules are disk shaped constituted of 7 (model 
2) tangent sphere of diameter a bonded together (figure 2). Model 2 is similar to the model 
considered by Edwards et. al. [^] with different kind of interactions. In case of model 1 the 7 
spheres forming the disc are hydrophobic and the 2 spheres at the two ends are hydrophilic. 



In case of model 2, the 6 outer spheres are hydrophihc and the inner sphere at the center of 
each molecule is treated as hydrophobic. Water molecule is modeled as a sphere. 

Attractive interaction between like particles (water-water, hydrophilic-water, and 
hydrophilic-hydrophilic) is modeled via a Lennard- Jones (LJ) potential: 

V{r)=Ae{{^r-{^f} (1) 

where r is the separating distance between the two interacting particles. The parameter 
e governs the strength of interaction and a defines a length scale. We have assumed that 
cr = 1 and e = 1 for all interactions. We have taken a cut-off r^ = 2.5cr, large enough to 
include excluded volume effects and attractive forces. 

Repulsive hydrophobic-hydrophilic and water-hydrophobic interactions are modeled by 
a truncated and shifted LJ potential (WCA potential): 

V{r)=Ae{(^f^-{^f} + e r <r,, 

= otherwise. (2) 

where r^ = 2^/^cr. 

For MC studies using LJ potential it is convenient to introduce reduced units. In terms of 
m, a and e, some of the quantities of interest are density p* = pa^, temperature T* = ksT/e, 
and pressure P* = Pa^/e. The simulation has been carried out both in NVT (fixed number 
of particles at constant volume and temperature) and NPT fixed number of particles at con- 
stant pressure and temperature) ensemble. During each MC step either chromonic or water 
molecules were chosen randomly and displaced using Metropolis criteria. The reorientation 
move was performed using quaternions 0. We have carried out simulations for different 
concentrations of the chromonic molecules at different temperatures. The concentration of 
chromonic is defined as 

where x is the number of atoms constituting a single chromonic molecules, N^ and Ny^ are 



the number of chromonic and water molecules respectively, x is 9 for model 1 and 7 for 
model 2. 

III. RESULTS 

In order to investigate the spontaneous formation of chromonic aggregates and their 
morphology, various choices of initial conditions were used. Model chromonic and water 
molecules are initially dispersed randomly in a L^ x Ly x Lz system with periodic boundary 
condition. Unless specified all of the results shown below are for model 1. Instantaneous 
snapshots of the columnar aggregates formed by chromonic molecules are shown in figure |^ 
and figure ^. At low concentration the chromonic molecules form short columns and with the 
increase in concentration the length and the number of the columnar aggregates increases. 
At higher concentration they form chain like aggregates. 

We have also calculated various positional and orientational pair correlation functions 
for the chromonic as well as water molecules to characterize the structure of the columnar 
aggregates, ^''^'^(r) is the radial distribution function of the center of mass of the chromonic 
molecules and is defined as 

^^^W = ;^(E^(r-r^l)) (4) 

g^^{r) is the same for solvent molecules (water) and is defined in an analogous way. gl^{r) 
and gf^ir) are the orientational correlation functions and are defined as 

9T{r) = ^^^^ (e A(ur ■ U-) 5(r - r^)^ , (5) 

9Tir) = -^^^ (e A(pur ■ pu-) Sir - r^)^ , (6) 

where r^| is the relative separation of the center of the ith and jth chromonic molecules, 
pu^*^ and u?'^ are a unit vectors lying parallel and perpendicular to the plane of the the ith 
chromonic molecule, P2{x) = 3/2cos^(x) — 1/2 is the second Legendre polynomial, and the 
sums range over all pairs of molecules. 



In practice, we work either with fully angle-averaged correlation functions, e.g. 

g''%r) = ^Jdng^'^{r,n), (7) 

or with cylindrically-averaged correlation functions, e.g. 

/'(r||,rx) = :^ rrf0/^(r||,r^,0). (8) 

zvr JO 

Here, ry and r± denote the components of r parallel and perpendicular to the orientation 
of the disc at the origin of the coordinate frame (see figure |lO|(a)), i.e. ry = r • u'^'^ and 
r_L = |r — (r ■ u'^'^)u'^'^|, and (p is the azimuthal orientation of r about u'^'^. g'^{r) , g^{r) and 
g'^{r\\,r_[_) are defined in an analogous way. 

In figure |^ we have plotted spherically averaged correlation function g'^'^{r) for different 
chromonic concentration and the first peak at r = 1.12cr signify strong positional correlation 
between chromonic molecules. Successive small peaks at r = 2.2(T, 3.3(J etc. indicate that 
the molecules aggregate in a columnar fashion as in figure |I|(a) or (b). 

The pair correlation between water molecules plotted in figure ^ exhibits the expected 
liquid like behavior with short range positional order. 

More detailed structural informations for the chromonic aggregate emerges from the two- 
dimensional correlation function 5''^^(r||, r^) and 5'i'^(^||, t±), shown in figure 0. With the refer- 
ence chromonic molecule situated at the origin the three prominent peaks at r^ ~ 1.1, 2.2, 3.3 
indicate that the chromonic molecules form columnar aggregates and the molecules within 
the columns is organized as in figure 0(a). Due to close packing effect the center of mass 
of the chromonic molecules within the columns can be shifted with respect to another by a 
which is reflected in the broadening of the peaks around ry ±1.0. Examination of g'^{r\\^ri_) 
reveals a strong intensity (~ 1.0) peak around ry ± 1.0 which indicates that the chromonic 
molecules within the columns are parallel to each other. 

We also compute average the aggregation number to study its dependence on chromonic 
concentration. We identify two chromonic molecules belonging to a cluster if the distance 
separating two monomers belonging to two chromonic molecules is < 2.5cr. From the cluster 
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size distribution we have calculated the average aggregation number through the following 
equation: 

<L>=Y.s'N{s)/Y.sN{s) (9) 

where N{s) is the number of cluster of size s. In figure § we have shown the average aggregate 
size < L > as a function of chromonic concentration. At low concentration < L >~ 0°'^. 
At higher concentration the dependence of < L > on is stronger than power law. Similar 
behavior is observed for conventional lyotropic systems [H,§. 

Since in our simulation we do not observe the formation of regular columnar aggregates as 
shown in figure |l|, we have simulated pre-assembled single and multiple columnar aggregates 
to test the stability of columnar aggregates. The simulation was done in NVT ensemble. 
The columnar aggregates remain stable over simulation time, in contrast to model 2 in 
which case pre-assembled columnar aggregates break and the system becomes isotropic. In 
addition our simulation with model 2 do not exhibits the spontaneous formation of columnar 
aggregates (figure |^). This indicates that we need to make judicious choice of the ratio of 
number of hydrophilic and hydrophobic units, an overall hydrophobicity of the chromonic 
molecules being necessary for the spontaneous self-assembly. 

IV. FREE ENERGY COMPUTATION BY UMBRELLA SAMPLING 

To calculate the free energy associated with the formation of columnar aggregates we 
have calculated the free energy of association in pulling apart two chromonic molecules by 
using umbrella sampling [|1^. The distance r^ separating two chromonic molecules were 



confined (figure |T0|) by means of a harmonic potential of the form : 

U{r) = ^k{r - r,f (10) 

where r^ is the equilibrium distance imposed between the center of mass of two chromonic 
molecules. The equilibrium distance were varied from l.lcr to la in steps of O.lcr - 0.2cr. 



Note that the equihbriuin distance between two chromonic molecules can not be less than 
a. From the individual biased average probabilities, P(r), obtained for each r^, unbiased 
free energy {F) can be constructed self-consistent ly using the weighted histogram analysis 



method (WHAM) |TT],|T2[. We have performed MC simulation at constant pressure for 25 
windows. The value of the spring constant used was k = 200. We have performed some 
simulation with higher values of spring constant k = 2000 and found qualitatively similar free 
energy behavior. The probability distribution of the distance r separating two chromonic 
molecules is , 

P(r) = {Sir - rj) 

= l|rfr^5[r-r,)]exp[-/3r(r^)], (11) 

for values of r for which P{r) is exceedingly small. Here, 5 is a Dirac delta function, N is the 
number of particles, r^ denotes the set of particle coordinates, /? = l/iksT), where ks is 
Boltzmann's constant and T is the absolute temperature, V^(r^) is the potential energy, and 
Z is a normalization factor (the configurational partition function). However, problem arises 
due to the fact that to P{r) becomes exceedingly small for values of r which give significant 
contribution to the free energy. Umbrella sampling makes use of a biasing potential to 
sample the region of phase space for which P[r) is exceedingly small. The distribution of r 
in the presence of a biasing potential is 

P\r) = ^^jdv''8[r-r,)] 

exp{-/3[r(r^) + f/(r)]} 
= |^exp[-/5f/(r)]P(r), (12) 

where Z' is the partition function for the biased Hamiltonian. From this it follows that 

Z' 
P(r) = — exp[/3f/(r)]P'(r). (13) 

Zj 

Thus, the distribution function Pir) can be obtained (to within a multiplicative constant) 
from a measurement of the biased distribution P'ir). The Helmholtz free energy F as a 
function of r can then be obtained (to within an additive constant) from 

8 



F(r) = -kBT\n[P{r)]. (14) 

By piecing together the relative free energies measured using a number of biasing potentials, 
it is possible to construct F{r) over any specified range of r. 

We have computed the potential of mean force of bringing two initially separated 
chromonic molecules on top of each other by the above method and the result is shown 
in figure 0. This case is referred to column 2 (figure |lO|(a)). From the free energy curve 
we see that there is strong attraction between chromonic molecules at short distance, lead- 
ing to face-to-face chromonic aggregation. Both energetic and entropic factors lead to this 
attraction: by stacking on top of each other the unfavorable hydrophobic interaction with 
water molecules is reduced and this stacking results in a net increase in the volume available 
to the water molecules and thereby increasing their entropy. At intermediate separation 
water molecules form solvation shells between the two chromonic molecules which results 
in a depletion repulsion between the chromonic molecules and thereby increasing the free 
energy barrier between them. At larger separation two chromonics do not interact with each 
other and so the potential of mean force levels off. However, at this stage we are unable to 
separate the entropic and energetic contribution in the free energy. 

We have also computed different orientational pair correlations gf^{r) and g2^{r) between 



the two chromonic molecules and are shown in figure |l2. 

From the gl^{r) we see that at the minimum separation (r = 1.0826cr) the two chromonic 
molecules are parallel and stacked together. However, their in-plane unit vectors make an 
angle to each other as is evident from g^{r). As the separation increases they start loosing 
their orientational correlations. 

Next we compute free energy of association of a chromonic molecules onto a colum- 
nar aggregates formed by two chromonics (referred to as column 3, figure |l^(b)) and three 
chromonic molecules (referred to as column 4, figure p^(c)) respectively. Within the colum- 
nar aggregate, the separation and orientation of the two chromonics are kept fixed which 
was computed at the location of the free energy minima for the column 2 case (see figure 



Til). The computed free energy is shown in figure 113. 



The free energy curves has the same features as in the case column 2 (figure |TT]). At large 
separation chromonic molecules do not see each other and consequently free energy reaches 
a plateau at larger separation. At intermediate separation the chromonic molecules start 
interacting with each other and start expelling the water molecules form the region between 
them to aggregate on top of each other and this leaves more room for the water molecules. 
This way water molecules gain entropy and chromonic molecules gain energetically (less con- 
tact with water molecules). This depletion effect eventually leads single chromonic molecules 
going on top of the existing columnar aggregates of length 2 or 3. Such type of depletion 
attraction and repulsion has been found between macromolecules immersed in a fluid of 



much smaller particles [|T^-|T5[. As in the case of column 2, we have also computed the 
orientational pair correlation for column 3 and 4. In case of column 3, the behavior is sim- 
ilar to the case of column 2. However for column 4, now the chromonic molecules remain 
perfectly parallel both within and perpendicular to the chromonic plane. The presence of 
three chromonic molecules within the aggregate strongly influence the fourth one whcih is 
coming on top of it. 

V. SUMMARY AND CONCLUSION 

We have developed a class of idealized model for the chromonic molecules. Simple site- 
interaction models represent a useful class of models for the study of the chromonic liquid 
crystalline self-assembly. For sufficiently hydrophobic molecular models formation of colum- 
nar aggregates is observed. Our free energy computation for different length of columnar 
aggregates gives valuable insight. However, more effort is required to test the isodesmic 
hypothesis. Future work will focus on improving configurational sampling (via rigid-body 
molecular dynamics) and on investigations of chemically realistic models of chromonic LCs. 
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FIGURES 



a 



fa) fb) fc) (d) (q) 

FIG. 1. Schematic representations of chromonic columnar aggregates, H-aggregates ( a and b), 
J-aggregates (c), brickwork structure (d) and layered brickwork (smectic) structure (e) 
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(2) 



Hydrophilic Unit (water like) 
( 1 ) ^ — - Hydrophobic Unit 

Water (Solvent) 

FIG. 2. A schematic representation of model chromonic and water molecules. 
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FIG. 3. Starting from an isotropic initial condition we get columnar aggregates. The parame- 
ters are: (f) = 0.081, Nc = 50, Nw = 5100, T* = 1.0 P* = 1.0, all in reduced units. For clarity the 
water molecules have not been shown. 
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FIG. 4. Same as in figure pi The parameters are: (j) = 0.137, Nc = 90, Nwr = 5100, T* = 1.0 
P* = 1.0, all in reduced units. For clarity the water molecules have not been shown. 
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FIG. 5. Pair correlation between the center of mass of the chromonic molecules for different 
chromonic concentration. 



ss . . 

g (r) 











^ = 0.0257 
^ = 0.109 




4 

r 



8 



FIG. 6. Pair correlation for water molecules for different chromonic concentration. 
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FIG. 7. Two dimensional correlation between the center of mass of the chromonic molecules 
(a) g'^'^{r\\,r±) and (b) gi^{r\\,r±). The parameters are same as in figure §. 
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FIG. 8. Average aggregation number < L > as a function of chromonic concentration 
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FIG. 9. Simulation with model 2 does not give rise to columnar aggregates. The chromonic 
density is (/> = 0.12. Other parameters are the same as in figure y. 
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(a) (b) (c) 

FIG. 10. Three different cases for wliich the free energy has been computed, (a) column 2, (b) 
column 3 and (c) column 4. 
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FIG. 11. Free energy as a function of the distance separating the center-of-mass of two 
chromonic molecules for the case UO^fa). 
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FIG. 12. Orientational pair correlations between two chromonic molecules for the case P^(a). 
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FIG. 13. Free energy as a function of the distance separating the center-of-mass of two 



chromonic molecules for the case shown in figure lC(b) and (c). 
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FIG. 14. Orientational pair correlations between chromonic molecules (a) for the case shown 



in figure 10(b) and (b) for the case shown in figure |l^(c). 
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